Quantum quasiresonances in grazing incident angle atom-surface collisions 
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The momentum transfer between the normal components to an index direction in the collision 
of an atom with a periodic surface is investigated. For fast atoms with grazing angle of incidence 
there is an interval of azimuthal angles around the index direction for which the energy transfer can 
be very efficient. This effect is reflected in quantum diffraction patterns with large non-specular 
peaks, associated with the parallel to the surface and normal to the index direction momentum 
component, and can be described in terms of quasiresonance. Although the classical dynamics 
does not reproduce the precise quantum diffraction probabilities, indicating the quantum nature of 
this effect, classical and quantum computations show that the span of the quasiresonance region 
coincides in both dynamics and can be classically estimated from the phase-space analysis. 

PACS numbers: 



Recent experiments in grazing angle of incidence col- 
lisions of fast atoms with a surface along a low index 
direction have shown quantum diffraction patterns even 
for a total translational energy as large as several keV 
[J El- El- This result is explained by the existence of 
two motion regimes, a fast one along a low index di- 
rection, with momentum p x — \J 2mE sin cos ip, and a 
slow one in the plane perpendicular to the index direc- 
tion, with momentum components p y = \/2mE sin 9 sin ip 
and p z = \J1mE cos 0, where E is the total incident en- 
ergy and m the atom mass. The component p z is chosen 
normal to the surface. For incidence along an index di- 
rection, ifi — and p y i — 0. Within these conditions, 
the fast atoms "feel" the surface as structureless parallel 
furrows aligned along the index direction [l[ , which con- 
stitutes an example of channelling on surfaces @, 0| ■ As 
a consequence the two regimes of motion are decoupled, 
resulting in a constant momentum component p x . Under 
grazing angle conditions the energy associated with the 
normal components is very small, and, as the change in 
p v is quantized, a quantum diffraction pattern is observed 
despite the relatively large total translational energy. 

However, the previous analysis does not explain to 
what extent the diffraction pattern still survives for inci- 
dence not in a surface index direction (p y i =/= 0), nor does 
it predict how efficient is the energy transfer between the 
normal components responsible for the diffraction pat- 
tern. A convenient analysis of the issue can be made in 
terms of quasiresonance. The concept of quasiresonance 
was first introduced to explain some very efficient and 
specific energy transfer between the vibrational and ro- 
tational degrees of freedom experimentally observed in 
state resolved collisions between vibrotationally excited 
diatom molecules and one atom 0| . An striking obser- 
vation was a remarkable insensitivity of the quasiresonant 
effect to the specific details of the interaction potential 



[S 0, HI • It was also shown that a fully classical treat- 
ment of the collision reproduces the experimental results 
remarkably well [1 0, Hi M, EH ■ Recently, the quasires- 
onance analysis has been extended to new processes in 
different classical systems II 2. If 3 1 , including grazing an- 
gle atom surface collisions [14j . 

An important open question is whether the quantum 
effects in the dynamics share any common mechanism 
with the classical quasiresonances. This issue has been 
analyzed in ultracold atom-diatom collisions [lo| . The 
results were not conclusive as for the low energies con- 
sidered the main inelastic quantum channels associated 
with quasiresonances remain closed. Grazing angle atom- 
surface collisions in a close-to an index incident direction 
provides a new ground to study the similarities between 
both dynamics. 

Classically, quasiresonance can be understood in terms 
of the adiabatic invariance theory and the method of 
averaging [la ]. Let us consider a two-dimensional inte- 
grable system described by the action variables (Ji, J2), 
perturbed by a interaction, depending on other addi- 
tional degrees of freedom, which is switched on and off. 
The interaction induces a coupling between the two "in- 
ternal" degrees of the perturbed system, and defines a 
process from an initial state (Ju,J2i) to a final state 
(Jif, ^2/)- Let us assume (a) a slowly varying perturba- 
tive interaction; and (b) an initial internal state satisfying 
the approximate M : N non-linear resonance condition 
Mu>2i — Nuju w 0, with M and N small integers, and 
u! a the two independent frequencies of the system. The 
analysis of dynamics of the nearly resonant system in 
terms of the standard secular perturbation theory intro- 
duces a canonical transformation to new action variables, 
h = -J\/N,h = -h + MJi/N, defined in a rotating 
frame in which the angle variable 02 associated with the 
action I2 oscillates much more rapidly than the other 
variables during the entire transient process. The aver- 
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aging of the transformed Hamiltonian over this rapidly 
oscillating angle leads to its conjugated action as the first 
term of the series expansion for the adiabatic invariant of 
the system in the proximity of the nonlinear resonance. 
In terms of the original actions, the existence of this local 
adiabatic invariance implies that 
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where A J — Jf — Ji. This accurate correlation between 
the action changes is the central result of the quasires- 
onance analysis [T3 , EH. The quasiresonance region is 
defined as the interval of initial values for which the con- 
ditions (a) and (b) are satisfied. 

In atom surface collisions, classical quasiresonance ap- 
plies to an effective two-dimensional free atom, described 
by the components of the momentum parallel to the sur- 
face, which is perturbed by the transient atom-surface 
interaction. Here the analysis can be simplified consider- 
ing a coordinate system with the x-axis fixed along the 
index direction of interest, since in this coordinate sys- 
tem the actions I\ and I2 can be directly identified with 
the components of the momentum p y and p x . The "ex- 
ternal" translational action that controls the strength of 
the interaction is p z . 

Assumption (a) corresponds to a collision with large 
translational energy (perturbative regime) and a small 
p z (slowly variant interaction); both together implying 
p z < < p x . Meanwhile the assumption (b) corresponds to 
p y « p x - These conditions define a scattering process 
with a grazing incident beam closely aligned to an index 
direction of the surface. Since x » y, z, the averag- 
ing over this rapidly oscillating variable (02 = x) allows 
us to identify its conjugated action variable, the momen- 
tum along such direction [1% = p x ), as the first order 
adiabatic invariant characterizing quasiresonance. The 
dynamics in the perpendicular plane {yz) is dictated by 
the Hamiltonian 



H m = 5- (pI+P 2 z ) +V av (y,z) 
2m 



with the averaged potential 



V av (y,z) = - 



V(x, y, z) dx, 



(2) 



(3) 



being L the surface periodicity along the x-axis. In a 
quasiresonance region, the coupling induced by this av- 
eraged potential enables an energy transfer between the 
normal components of the momentum, p y and p z - The 
efficiency of such energy transfer is determined by the 
form of V av resulting from the original potential and the 
particular index direction implied. In general, the energy 
transfer is more efficient for incidence close to a low index 
direction. Outside quasiresonance such energy transfer 
does not occurs, and incident beams "feel" a plain sur- 
face and only suffer an effective specular reflection. 



In a coordinate system with the x-axis fixed along the 
surface direction (1,0), the adiabatic invariance of the 
momentum along a symmetry direction (M, N) , associ- 
ated with the approximated resonance condition Mp' y — 
NjpL ~ 0, leads to the correlation Ap' y /Ap' x = -M/N 
14j . In this coordinate system the efficient energy trans- 
fer in the plane perpendicular to the direction (M, N) is 
reflected by significant changes in both internal actions, 
p' x a.nd p' y . This is the usual result of the general quasires- 
onance analysis. 

Quasiresonance is satisfied for grazing incident beams 
aligned over an interval of azimuthal angles around the 
crystal index direction. We define the width W py of the 
quasiresonance region as the interval of initial values of 
p y for which the approximated resonance condition (b) is 
satisfied, or, alternatively, W v as the corresponding in- 
terval of incident azimuthal angles tfi. Both quantities 
are related by W v = 2 sin -1 [W py /(2\ / 2mE sin 6i)\ . W py 
can be estimated from the analysis of the phase space as- 
sociated with the averaged Hamiltonian at frozen values 
z f r of the normal z coordinate [HI • We focus on the dy- 
namics around the resonant incident direction, p y i = 0, 
which can be properly described by the pendulum like 
Hamiltonian H nPn d = P y /2m + V av (y, Zf r ). The maxi- 
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mum excursion in action of the islands of the resonance 
zone in this integrable system provides the simple expres- 
sion 



W py w 2^2m [V av (y u 



z m ) - V av (y s ,z m )], (4) 



being z m the normal coordinate that gives the maximum 
strength interaction at incident energy E, y u the hyper- 
bolic point associated with the innermost separatrix, and 
y s the most stable fixed point of the potential V av (y, z m ). 
The expression Q is valid provided that the resonance 
zones associated with the different crystal index direc- 
tions remain well isolated in phase space. In a non per- 
turbative regime, the overlapping between neighbouring 
resonance zones in phase space makes the quasiresonance 
regions shrink [l3j ]. 

To illustrate quasiresonance we simulate the atom- 
surface collision considering a Lennard- Jones and Devon- 
shire potential with parameters taken from Ref. [lij ]. 
describing the elastic scattering of He atoms from a per- 
fectly periodic rigid LiF(OOl) surface. The low normal 
velocities implicit under grazing angle conditions, and 
the robustness of the quasiresonance phenomena with re- 
spect to the details of the interaction potential justify the 
use of this simple model. 

Figure [T] shows the changes in p x and p y for classical 
trajectories closely aligned to an index direction. Two 
different behaviors are clearly distinguished; there are 
significant changes in p y in the region around the sym- 
metry direction <pi = 0, identified as the quasiresonance 
region, and a constant value of p y outside it. The adia- 
batic invariant in the quasiresonance p x remains constant 
in the whole depicted range. The change in p z (not shown 
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FIG. 1: (Color online) Ak x = (p x f —p X i)/h (black dashed- 
dotted line) and Ak y (brown dots) versus the initial azimuthal 
angle (pi, for classical trajectories closely aligned to the [110] 
direction (ipi = 0). All the trajectories have incident polar 
angle 6i = 0.5067T, initial Zi = 60A and (xt,yi) chosen at 
random in a unit cell. The red solid line and the blue dashed 
line represent the mean values and the root-mean-square val- 
ues of Aky respectively, obtained at each angle tfii from an 
average over 2000 classical trajectories with identical initial 
conditions, except for the coordinates (xi,yi). The mean val- 
ues (red squares) and the root-mean-square values (blue cir- 
cles) derived from the quantum diffraction probabilities are 
also depicted. 



in the figure) is dictated by the energy conservation in 
the global elastic process. 

For the quasiresonance conditions (a) and (b), the 
quantum mechanical treatment of the collision can be 
simplified, as only a small number of quantum channels 
are effectively open [l| . Under these conditions the large 
momentum component along the index direction remains 
constant and only quantum transitions corresponding to 
reciprocal lattice vectors perpendicular to this direction 
present a non-negligible diffraction probability. As only a 
small amount of energy associated with the perpendicular 
degrees of freedom is available, only a few quantum chan- 
nels contribute to the diffraction spectrum. This effect 
has been experimentally confirmed for incidence along a 
low crystal index [H, Then the quantum dynamics 

can be described by an effective 2D model [lj with the 
averaged potential Eq. 

Figure [5] shows the diffraction spectra for four inci- 
dent azimuthal angles: along a crystal index direction, 
inside, on the edge, and outside the classical quasires- 
onance region. The significant non-specular diffraction 
peaks in the first three spectra correspond with the effi- 
cient energy transfer between the perpendicular momen- 
tum components, p y and p z , predicted in the classical 
quasiresonance. As the incident direction approaches the 
edge the intensity of the inelastic channels decreases, be- 
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FIG. 2: Quantum (solid lines) and quasiclassical (dashed 
lines) diffraction spectra for a grazing incident beam with 
energy E — 200 eV, polar angle 9i = 0.5067T and different 
azimuthal angles (a) along the [110] direction, ipi — 0, 
(b) inside the quasiresonance region, tfii = 0.01, (c) on the 
edge of the quasiresonance region, ipi = 0.02, and (d) outside 
of the quasiresonance region, ipi — 0.04. The quasiclassical 
probabilities were calculated replacing the continuous distri- 
bution by an histogram. The diffraction probabilities have 
been convoluted with a Lorentzian function to account for 
possible broadening effects. 



ing the specular diffraction peak the only that survives 
for incidence outside the quasiresonance region. The pre- 
cise distribution of the diffraction probabilities is dictated 
by quantum effects in the dynamics of the system, and 
as expected, is not correctly reproduced by the classical 
model, see figure [2j The important result however is the 
qualitative agreement, and especially the correct predic- 
tion of a significant qualitative change in the quantum 
diffraction spectrum as the incident direction crosses the 
border of quasiresonance. Both the exact classical and 
exact quantum results confirm the simple quasiresonance 
analysis. As figure [T] shows, there is also a very good 
agreement between the trends in classical and quantum 
changes in p y . Hence both classical and quantum dynam- 
ics reflect the same mechanism of efficient energy transfer 
that characterizes quasiresonance around the low index 
direction. 

Figure [3] shows the dependence of the span of the 
quasiresonance region on the total translational energy. 
The monotonic increase of the width W Py with the energy 
is a peculiar feature of the simple interaction potential 
considered. Since in this model the averaged potential is 
zero except for the principal surface directions (1,0) and 
(0, 1), quasiresonance around other symmetry directions 
is a higher order perturbative effect. Therefore, the main 
mechanism for the eventual demise of the main quasires- 
onance regions, i.e., the overlapping between neighbour- 
ing resonance zones in phase space, is nearly absent. The 
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FIG. 3: The span of the quasiresonance region around the 
[110] direction versus the incident energy E, in a process with 
grazing incident angle 8i — 0.5067T. Labels on the left corre- 
spond the angular width W v ; the squares are the numerical 
classical estimation, the crosses the numerical quantum es- 
timation, and the solid line the classical prediction derived 
from 0. Labels on the right give the size of the quasireso- 
nance region in momentum space, Wk v = W py /h; the dashed 
line corresponds to the classical estimation ((4]). As a cri- 
terion for the numerical estimation of the angular width of 
the quasiresonance region, we consider the full width at half 
height of root-mean-square values of Ak y . 



angular with W v decreases rapidly at low energies and 
changes very slowly at higher energies. At low energies, 
no normal quantum inelastic channels are energetically 
allowed and, as expected, there is not good agreement 
between the quantum and classical dynamics. Both re- 
sults become nearly identical as the energy increases and 
more diffraction channels are open. 

In conclusion, the grazing angle collision of an atom 
with a periodic surface for incidence closely aligned to a 
(low) index direction has been investigated using both a 
classical and a quantum description. For a large interval 
of incident energies the results can be analyzed in terms 
of quasiresonance. This interval is limited below by low 
energies for which the dynamics is not perturbative, and 
above by large energies with large normal to the surface 
component of the momentum for which the atom-surface 
interaction is not slowly variant. Quasiresonance is char- 
acterized in this system by the adiabatic invariance of the 
momentum component along a symmetry direction of the 
target surface, together with an efficient energy transfer 
between the normal degrees of freedom. Even when the 
classical results do not reproduce the precise quantum 
diffraction probabilities, the averaged results are in good 
agreement in both dynamics, specially for large energies. 
The width of the quasiresonance region can be classically 
estimated from the evolution of the phase-space of the 
perturbed system in the transient process. The angular 
width W v decreases rapidly at low energies and remains 
nearly constant for large energies, provided that the in- 
teraction does not activate high order quasiresonances. 



We have also analyzed the analogies between the clas- 
sical and quantum quasiresonance behavior. The mayor 
disagreement between both dynamics occurs at low ener- 
gies, due to the energy threshold of the quantum diffrac- 
tion channels. As was described in ultracold atom- 
diatom collisions [Io| . no quantitative similarity is found 
when the inelastic channels are closed. For large trans- 
lational energies, the quasiresonance width coincides, in- 
dicating a common mechanism in both dynamics. 

Our results for energies up to hundreds of eV indicate 
an angular quasiresonance width W v « 0.02 rad, which is 
above the beam divergence achieved in experimental se- 
tups. Although we have assumed a perfect rigid surface, 
this is a good approximation for large beam energies as 
the particle wavelength is much smaller that the thermal 
vibration of the lattice atoms [l|. Therefore, experimen- 
tal investigations to characterize quasiresonance regions 
in more realistic atom-surface interactions, measuring the 
diffraction pattern for incident beams closely aligned to 
the low crystal index direction are in principle possible. 
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